Open AccessCCS ChemistryRESEARCH ARTICLES22 Oct 2022Structure Identification for Force-Induced Reaction Using Single-Molecule Conductance Measurement Jueting Zheng†, Wenli Gao†, Taige Lu, Lijue Chen, Luchun Lin, Ruiyun Huang, Yongxiang Tang, Gang Dong, Junyang Liu, Yifei Pan, Wengui Weng and Wenjing Hong Jueting Zheng† State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering and Institute of Artificial Intelligence and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen 361005 †J. Zheng and W. Gao contributed equally to this work.Google Scholar More articles by this author , Wenli Gao† Department of Chemistry, College of Chemistry and Engineering, Xiamen University, Xiamen 361005 †J. Zheng and W. Gao contributed equally to this work.Google Scholar More articles by this author , Taige Lu State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering and Institute of Artificial Intelligence and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen 361005 Google Scholar More articles by this author , Lijue Chen State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering and Institute of Artificial Intelligence and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen 361005 Google Scholar More articles by this author , Luchun Lin State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering and Institute of Artificial Intelligence and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen 361005 Google Scholar More articles by this author , Ruiyun Huang State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering and Institute of Artificial Intelligence and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen 361005 Google Scholar More articles by this author , Yongxiang Tang State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering and Institute of Artificial Intelligence and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen 361005 Google Scholar More articles by this author , Gang Dong State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering and Institute of Artificial Intelligence and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen 361005 Google Scholar More articles by this author , Junyang Liu State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering and Institute of Artificial Intelligence and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen 361005 Google Scholar More articles by this author , Yifei Pan Department of Chemistry, College of Chemistry and Engineering, Xiamen University, Xiamen 361005 Google Scholar More articles by this author , Wengui Weng *Corresponding authors: E-mail Address: [email protected] E-mail Address: [email protected] Department of Chemistry, College of Chemistry and Engineering, Xiamen University, Xiamen 361005 Google Scholar More articles by this author and Wenjing Hong *Corresponding authors: E-mail Address: [email protected] E-mail Address: [email protected] State Key Laboratory of Physical Chemistry of Solid Surfaces, College of Chemistry and Chemical Engineering and Institute of Artificial Intelligence and Innovation Laboratory for Sciences and Technologies of Energy Materials of Fujian Province (IKKEM), Xiamen University, Xiamen 361005 Google Scholar More articles by this author https://doi.org/10.31635/ccschem.022.202202297 SectionsSupplemental MaterialAboutAbstractPDF ToolsAdd to favoritesDownload CitationsTrack Citations ShareFacebookTwitterLinked InEmail Spiropyran derivatives are prototype mechanophores with a promising application as molecular sensors because of their changeable structure under external force stimuli. However, the chemical structure evolution under external stimuli remains unclear due to the uncertainty and difficulty in distinguishing the structures of different ring-opened merocyanine isomers generated in the force-induced reaction. Here we identify the structure of isomers produced by the force-induced reaction of spiropyran derivatives using a single-molecule conductance measurement and an unsupervised clustering algorithm. We found that the original data from the single-molecule conductance measurement can be divided into four clusters through unsupervised clustering. By introducing a photoinduced reaction and theoretical calculation, we identified and attributed the four clusters of data to the multiple states of the molecular junctions. Our work demonstrates that a single-molecule break junction measurement can distinguish the isomers in the force-induced reaction, suggesting the great potential of single-molecule conductance measurement and unsupervised clustering approaches for structural analysis. Download figure Download PowerPoint Introduction Mechanochemistry offers a unique strategy for the study of mechanical transformations of bonds in mechanophores.1–6 Among these mechanophores, spiropyran-cored (SP) compounds can not only be mechanically activated but also stimulated by other cues to display color changes due to the spiro C–O bond breaking into merocyanine (MC) isomers. Therefore, spiropyran derivatives have been extensively studied, especially for the incorporation into polymers and the corresponding polymer mechanochemistry.7–10 Investigation of the evolution of MC isomers under external stimuli is critical to understanding the structure–property relationships of MC and the application of SP mechanophores.11–13 However, since there are several product isomers generated in the force-induced reaction, the exact chemical structure of spiropyran derivatives formed under external stimuli remains unclear.14 The in situ characterization of the structural changes in spiropyran derivatives under external stimuli has the potential to reveal the structural evolution in the forced-induced reaction process.15 Single-molecule techniques16–18 hold great potential to reveal the structural evolution of mechanophores in the force-induced reaction process. Atomic force microscopy (AFM)-based single-molecule force spectroscopy has been adopted to study SP-containing polymers.19,20 However, the structure of the MC isomers13 generated in the force-induced reaction cannot be well distinguished due to the presence of multiple isomers in the polymers. Single-molecule break junction techniques, including mechanically controllable break junction,21–23 scanning tunneling microscopy break junction (STM-BJ),16,17,24–28 and atomic force microscopy break junction,29,30 can also analyze the structure evolution in a reaction at the single-molecule level,31–33 which offers a unique chance for investigating the structural changes during the reaction of mechanically sensitive compounds by monitoring the single-molecule conductance changes under applied tension force. In this study, we synthesized a novel thioether-functionalized spiropyran-cored compound ( SSP) and measured its single-molecule conductance using the STM-BJ technique with repetitive force applied to the junction. We divided over 4000 individual curves into four types via a lab-developed unsupervised spectral clustering algorithm according to the Calinski–Harabasz (CH) index.34–37 The change of conductance along with the stretching distance indicates that the Au/ SSP/Au junction experienced structural changes during the break junction processes. STM-BJ measurements on the merocyanine compound ( SMC-UV), formed from ultraviolet irradiation of SSP,38 and a reference molecule SIN-Ref further suggest that the highest conductance value among the four types originated from a shorter configuration. We attributed the other three types to the ring-closed form of SSP and possible merocyanine ( SMC) products based on the experimental observations and clustering algorithm combined with density functional theory (DFT) calculations. Our results demonstrate the potential of single-molecule break junction measurement and clustering algorithm for structural analysis of the chemical reactions. Experimental Methods Materials synthesis 4-Iodophenylhydrazinium hydrochloride, iodomethane, 6-ethynyl-4,4-dimethyl 2,3-dihydrothiochromene, 2-hydroxy-3-nitrobenzaldehyde, PdCl2(PPh3)2, and CuI were purchased from Aldrich (Shanghai, China). 3-Iodo-4-nitrophenol and methyl isopropyl ketone were purchased from Meryer (Shanghai, China). All the other chemical reagents were purchased from Sinopharm (Shanghai, China) and used without further purification. The detailed synthesis procedures and characterization data (1H NMR, 13C NMR, and mass spectrometry) of SSP and SIN-Ref compounds are listed in the Supporting Information Figures S1–S3. Data analysis details Conductance histograms were constructed from the conductance traces with a bin size of 1000 (bin width = 0.01 log(G/G0)). Two-dimensional conductance-distance histograms were constructed with a bin size of 1000 × 500 (log(G/G0) × nm). The conductance-distance traces were aligned by taking Δz = 0 for the conductance value equal to 10−0.3G0. The spectral clustering algorithm was adopted from the previous report.37 The data analysis was processed in our open-source code pace XME-data analysis (Xiamen Electronics, https://github.com/Pilab-XMU/XMe_DataAnalysis). Here, each conductance trace was constructed into the conductance histogram. By choosing conductance as the feature, spectral clustering with a CH index was applied to assign similar conductance histograms, thus dividing the original conductance traces into clusters. Theoretical calculation details DFT calculation was performed with the B3LYP/6-31G* in Gaussian 16 package.39 The transmission spectra were performed using the QuantumATK 2018.6 package based on the DFT-nonequilibrium Green’s function (NEGF) method.40 Results and Discussion Single-molecule conductance measurements We fabricated single-molecule junctions via a lab-built STM-BJ setup,41 as shown in Figure 1a (for experimental details, please refer to Supporting Information Section 2 and Figures S14–S15). Since the force between the Au and S atom can transmit to the core of the target structure, the break junction process formed a ring-opened isomer of SSP compound ( SMC). This mechanochemistry reaction has several possible SMC isomer products (Figure 1b),13 thus the measuring results of SSP junction might contain multiple types of conductance traces. To classify the different states of the junction, we utilized an unsupervised spectral clustering algorithm to separate all curves of Au/ SSP/Au junction tests.37 As shown in Figure 1c, the original conductance histogram did not show a noticeable conductance feature, and the conductance traces exhibited a variety of the plateau position before clustering (inset). After analyzing the original data with a spectral clustering algorithm, the CH index42 presented the highest value at four clusters, indicating that the data is most likely separated into four types, as demonstrated in Figure 1d. The original data (4580 traces) can be divided into four kinds of typical conductance traces with distinct plateau features at certain conductance ranges after clustering, as shown in Figure 1e inset. The clustered traces were statistically analyzed into conductance histograms, which are illustrated in Figure 1e. The trace amounts from clusters 1, 2, 3, and 4 are 868, 1073, 1209, and 1430, respectively. Figure 1 | (a) Schematic of STM-BJ technique. (b) Left panel: the chemical structures of SSP and its ring-opened form SMC. Right panel: possible structures of the isomers of SMC. (c) The conductance histogram of the SSP junction without data selection. Inset: typical conductance traces before clustering. (d) Calinski–Harabasz (CH) index of the clustering results for the SSP junction and schematic of spectral clustering. (e) The conductance histograms constructed from the traces according to the clustering algorithm. Inset: clustered conductance traces. Download figure Download PowerPoint 2D conductance-distance analysis After clustering, we found an abnormal “jumping” feature in many of the conductance traces in clusters 2–4, as shown in Figure 2a. That is, when the distance of electrodes increases, the conductance value of the junction is enhanced. However, the transmission probability should be lower when the electrode distances are larger. Walkey et al.15 observed the in situ junction conductance change for spiropyran derivatives with Au–C bond connection in the individual conductance-distance traces, whereas the direct cause of the in situ changed conductance state remains unknown. Our results also indicate a structural change to a more conductive form during the break junction process. To further investigate the correlation between conductance and the structure of single-molecule junctions, we applied a 2D conductance-distance histogram43 to obtain the conductance evolution during the stretching process. As shown in Figure 2b, we selected those curves with this anomalous “jumping” feature and constructed them into another conductance-distance histogram containing 995 curves (out of 4580 curves).44 This phenomenon indicated that the structures of the single-molecule SSP junctions changed during the stretching process, leading to the conductance increase rather than decrease over a long distance. To investigate the new structure generated in the stretching process, we constructed the 2D conductance-distance histograms of the four clusters in Figure 2c–f. The corresponding 2D conductance-distance histograms constructed from the curves with the “jumping” feature are shown in Supporting Information Figure S16. Clear plateaus demonstrate the formation of molecular junctions with various conductance features and their corresponding distances. Among these four clusters, Cluster 1 has the shortest plateau, indicating that the junction length of cluster 1 is smaller than the other three. In all, 2D conductance-distance histogram analysis further confirmed the structural changes in the SSP junction, and the break junction process contained multiple configurations with different lengths. Figure 2 | (a) Typical conductance-distance traces with “jumping” feature. These traces were treated with offsets in the x direction for better view. (b) 2D conductance-distance histogram constructed from the curves with conductance enhancement with the stretching. The traces were realigned at the point where the conductance value changed abruptly. For the realignment, we chose the middle of the x-axis scale (Δz = 1.5 nm) for better vision. (c–f) 2D conductance-distance histograms constructed from the four clusters. Insets: length distributions of the molecular junction for the corresponding clusters. Download figure Download PowerPoint Structure identification To reveal the structure of the shorter junction configuration in the measurements of SSP junctions, we applied 365 nm UV-light irradiation to form the ring-opened merocyanine compound SMC-UV ( Supporting Information Figures S17 and S18). We also designed and synthesized the reference molecule, SIN-Ref, which contains the indoline structure but is half the length of the SSP and SMC molecules, as shown in Figure 3a. We measured SMC-UV and SIN-Ref junctions with the same methodology and observed an apparent plateau in conductance traces. There is also a multi-conductance feature for the SMC-UV junctions, shown as the purple line in Figure 3b. The high conductance is located around 10−3G0, and the low conductance is around 10−4.7G0. The high conductance value of spiropyran derivatives has been reported and attributed to the force or UV light-activated ring-opened configurations.15,45 However, such attribution cannot explain the origin of the low conductance feature in our experiments. The high conductance feature of the SMC-UV junction is similar to that of Cluster 1 in SSP and SIN-Ref. The 2D conductance-distance histograms of SMC-UV (Figure 3c) and SIN-Ref (Figure 3d) show similar features in the region around 10−3G0. The corresponding distance distribution (insets) demonstrate that the junction length is consistent, indicating a similar configuration in three individual measurements. The broad conductance peak constructed from the SSP-Cluster 1 indicates that a configuration of this type of junction is unstable, resulting in the decentralized distribution of conductance.32 The comparison demonstrates that the higher conductance value for the SSP and SMC-UV junction originates from the junctions with a similar configuration to SIN-Ref rather than the fully-stretched, higher conductance merocyanine form. Figure 3 | (a) The chemical structures of SMC-UV and SIN-Ref. (b) Conductance histograms of SSP-Cluster 1, SMC-UV, and SIN-Ref. (c) 2D conductance-distance histograms of SMC-UV. Insets: the corresponding distance distributions. (d) 2D conductance-distance histograms of SIN-Ref. Insets: the corresponding distance distributions. (e) The configuration demonstrations of high and low conductance. (f) Conductance histograms of the clustered results in the low conductance region of SSP and SMC-UV junctions. Download figure Download PowerPoint According to the above analysis, the junction formation configuration of the SSP and SMC-UV junction can be divided into two processes: First, when the gold electrodes separate and the distance increases, half of the SSP/SMC-UV molecule is captured through a Au–N connection,46 resulting in a high conductance with similar short plateau length distributions to the SIN-Ref junction; second, the distance between the gold electrodes increases gradually, the short junction is ruptured, and the complete molecule can be determined in the low conductance region. The possible molecular configuration evolution is shown in Figure 3e. Besides Cluster 1, there are three other types of curves extracted from the original data of the SSP junction in the low conductance region of 10−3.5 to 10−5.0G0 from the clustering, suggesting that the different structures in the isomers could be distinguished using the single-molecule break junction techniques. To further understand the origin of the conductance values, we performed spectral clustering of SMC-UV and compared the results with SSP, as shown in Figure 3f. The multi-conductance feature from the clustered SMC-UV junction agrees with the SSP junction in the low-conductance region, indicating similar structures of the molecules in these two junctions. We found that the inapparent lower conductance peak of the SMC-UV junction was composed of three different types of conductance traces. This phenomenon was probably caused by the part of the SSP molecule that did not transform into the merocyanine form after UV irradiation due to the spiropyran-merocyanine equilibrium.38,47 By comparison, molecular junction configurations with notably different lengths in the SSP junctions were distinguished, whereas the origin of the three clusters in the lower conductance regions with similar lengths remained unclear. Because the energy barriers for the conversion between the possible isomers of spiropyran derivatives are low,48,49 complex products appear in the ring-opening processes. Therefore, there are multiple compounds in the break junction processes of both SSP and SMC-UV junctions so that the molecular structure of the configurations needs further identification. DFT-NEGF calculations To reveal the possible structures corresponding to the different conductance values, we conducted theoretical calculation to get transmission spectra of the SSP compound and its four SMC isomers of ring-opening reaction, TTT, TTC, CTT, CTC, using DFT combined with the NEGF.40 According to the previous reports, the product SMC has several isomers, as shown in Figure 4a.13,50 Our calculation revealed that the ring-closed SSP compound is less planar than the ring-opened SMC, as shown in Figure 4b ( Supporting Information Figure S19), which also agrees with the previous reports.51,52 The conjugation was enhanced as the planarity changed from SSP to SMC, indicating that the ring-closed SSP junction is less conductive than SMC. As shown in Figure 4c, the transmission order of these isomers from low to high is SSP< CTC≈ CTT< TTC≈ TTT. From the transmission spectra, the ring-closed SSP has the lowest conductance among the different structures, and when it becomes more conjugated from ring-opening, the charge transport will be enhanced.53 This is also the origin of the “jumping” feature for the SSP junction shown in Figure 2b. Figure 4 | (a) The chemical structures calculated by DFT-NEGF. (b) The configurations of the SSP, CTC, and TTC isomer junctions (c) Transmission spectra of the closed and the ring-opened isomers. Download figure Download PowerPoint Combining the results of single-molecule conductance and DFT calculation, we can infer that the conductance values of the SSP junctions, which locate at 10−4.5 and 10−4.2G0, are originated from the ring-opened isomers. We attribute the three types of conductance traces with the low-conductance value from the SSP junctions into TTC/ TTT (cluster 2), CTC/ CTT (cluster 3), and closed-form SSP (cluster 4), from high to low, according to the results of conductance measurement and DFT calculations. According to the small energy barriers of isomerization of CTC/CTT and TTC/ TTT,13,48 we could not further classify the two structures in the two groups due to the interconversion between the two structures at room temperature. The proportion of isomers can also be estimated from the calculated free energy of each isomer. TTC was reported to be the most stable configuration among MC forms under thermal54 and UV-light55 stimuli, indicating that it should be the most common form. Interestingly, in our results, the number of curves from TTC/TTT junctions (cluster 2 containing 1073 curves) is close to that of CTC/CTT junctions (cluster 3 containing 1209 curves), suggesting that the structure distribution of the different isomers in a force-induced reaction could change from the equilibrium state in the confined nanogap of the single-molecule conductance measurement. Conclusion We investigated the in situ force-induced structure transition for thioether-functionalized spiropyran ( SSP) derivative at the single-molecule scale via charge-transport characterization using the STM-BJ technique. We utilized an unsupervised clustering algorithm to classify the several types of traces in the complicated processes during stretching and reforming the SSP junction. Among the clustered four types of conductance traces according to the CH index, one cluster displayed a short plateau feature, whereas the other three presented longer plateaus. By introducing UV irradiated form ( SMC-UV) and a reference molecule SIN-Ref, we revealed that the cluster with a short plateau with the highest conductance originated from an incompletely stretched configuration. In addition, the other clusters with longer plateaus can be attributed to SSP and the ring-opened isomers compound ( SMC) according to DFT-NEGF calculations. Our results revealed the origin of different conductance features at the single-molecule scale, suggesting that the structure of different isomers generated in the force-induced reaction can be extracted, demonstrating the great potential of the single-molecule junction technique in the structural analysis of chemical reactions. Supporting Information Supporting Information is available and includes synthesis and characterization, conductance measurement details, and theoretical calculation details. Conflict of Interest There is no conflict of interest to report. Funding Information This work was supported by the National Natural Science Foundation of China (grant nos. 22173075, 21933012, 61901402, 31871877, and 21774106), the National Key R&D Program of China (grant no. 2017YFA0204902), the Fundamental Research Funds for the Central Universities (grant nos. 20720200068 and 20720190002), and the Natural Science Foundation of Fujian Province (grant no. 2018J06004).